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ABSTRACT 

Fusion reaction raies for DT, DD and DHe^ reactions 
have been calculated for non-Maxwe Ilian distribution of 
velocities. In particular, <<3v> values at various kinetic 
temperatures have been calculated for a two - tei^erature 
pseudo — Maxwellian distribution, which may approximately 
represent the velocities in a magnetic piirror reactor. 

It is found that for the same kinetic temperature, 
a two - temperature distribution gives enhanced reaction 
rates betv/een 5 to 15 keV for DT reaction and between 20 to 
40 keV for DD and DHe^ reactions. These are the tempera t-ure 
ranges of maximum practical interest. Calculations at other 
temperatures are also reported. It is found that at very 
high temperatures, the two - temperature distribution gives 
reduced reaction rates in contrast to the favourable beha- 
vioiar at lower temperatures. 
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ft 
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^1 ) 
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Noimaj-ized ifexv/e Ilian distribution function 
of velocity for species A 

%( ^2 ) Normalized Maxwellian distribution fionction 
of velocity for species B 
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Chapter 1 


INTRODUCTION 

As is v/ell knovm, in both the lightest and the 
heaviest nuclei, the mean binding energy per nucleon is less 
than that in the intermediate ones. Consequently, nuclear 
energy is released in the fission of heavy elements and in 
the fusion reactions of the lightest ones. The fusion x’eac- 
tions of practical interest at present in the nuc 3-ear fusion 
program are the following ; 


D + T ^ He^"^ (3.5 MeV) + n (14.1 MeV) (1.1) 

D + D ^ T (1 MeV) + p (3 MeV) (l.aa) 

D + D > He^(0.8 MeV) + n (2.5 MeV) (1 .2b) 

5 L 

D + He ^ He^ (3.7 MeV) + p (14.6 MeV) (1.3) 


For the above reactions to occur, the two reactant 
particles must overcome the Coulomb barrier. In thermo- 
nuclear fusion this is achieved by raising the temperature 
of the reacting species. For useful reaction rates, the 
temperatures required are 5 to 10 keV for the DT reaction, 
and 20 to 40 keV for the DD and DHe^ reactions [l]. These 

* O 

are extremely high temperatures, of the order of 10 K, and 
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most of the difficulties of achieving nuclear fusion are 
related to this. In contrast, we note that a fission reac- 
tion hetween a neutron and, say, a uranium nucleus, can 
proceed at room temperature, because neutron is not a charged 
particle. 

The reacting mixtures in Eqs. 1.1 to 1.3 become fully 

Q 

ionized plasma at temperatures of the order of 10 K, Using 
the cross sections for the reactions under consideration, the 
reaction rates can be calculated if the velocity distributions 
of the reacting species are known. The reaction rates are 
generally calculated by assuming Maxwellian distribution of 
velocities of the two reacting species [lj. 

Our purpose in this work is to calculate the reaction 
rates of the reactions listed in Sqs. 1.1 to 1.3 if the 
velocity distributions deviate from the Maxwellian. This can 
be considered of some practical interest in view of the fact 
that velocity distributions in fusion reactors may be generally 
considered as different from Maxwellian, even though close to 
it. 

For this purpose, we analyze in detail the methods of 
reaction rate calculations in Chapter 2, and develop appropriate 
equations for calculating the reaction rates under no n-M>Iaxwe Ilian 
velocity distributions. Detailed calculations are carried out 
for a two - temperature pseudo - Maxv/ellian distribution. The 
results are reported and discussed in Chapter 3. 



Chapter 2 


ANALYSIS FOR REACTION RATE CALCULATIONS 


2.1 THE BASIC EQUATIONS FOR BINARY REACTION RATES 


We consider a system containing any nimiber of species 
of particles, and direct our attention on tvjo particular 
species A and B, and infact on a particular binary reaction 
between these two species. Let n^ and n^ be the number 
densities per unit volume of the two species, respectively, 
at a particular point under consideration, and f y., ( v^ ) and 
^B ^ ^ their normalized velocity distribution functions 

at the same point, i.e., 


A 


f: 








1 ( 2 . 1 ) 


3 3 

Here d v^ and d V2 are infinitesimal volume elem.ents in the 

velocity spaces. For examxjle, in cartesian coordinates, 

d^v^ = dv^^ *^^y1 ' '^"'^zl * Choosing v^ for the velocity space 

> 

of species A, and v^ for that of species B is for the conve- 
nience of subsequent discussion, in which one wants to 
distinguish between the velocities of two particles even if 
the species A and B are taken to be the same. 

Let a(v) be the microscopic cross section for the 
binary reaction mder consideration, lArhere 
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V 



( 2 . 2 ) 


If one considers, as a target, a single particle of species 
A moving with velocity v^ , and n^fg (v^) d'^V 2 | V 2 -v^ | = 

d^Vg.v as an incident beam of species B, the result- 
ing removal rate from the beam v/ill be 


( V2 ) d"v2 V a(v) . 


(2.3) 


^ 'y 

Multiplying Eq, 2.3 by n^f^ (v^) d v^ , the number of particles 

) 

of species A per unit volume with velocities betv/een v^ and 

— ^ ^ / 

v^ + dv^ , and integrating over the respective velocity spdces, 

we obtain 


R = 


^A^B 


J 



d^Vg V o 


(v) 




"B 




(2.4) 


where R is the number of interactions per unit volume per unit 
time of the desired type, and v is as in Eq. 2.2 


An alternative way of obtaining (or looking at) Eq, 2.4 

f — 3 

is to consider (v^^ d-^v^.v as the relative path length 

travelled by particles of species A (between v^ and v^ + dv^) 
in a unit spatial volume per unit time, and to consider 
a(v) Hgfg ("^2^ d^V 2 as the macroscopic cross section due to 
particles of species B betv/een and v^, + dv 2 . Since the 
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latter gives the probability of interaction per unit (relative) 
path length travelled, multiplying the two and integrating over 
the velocity spaces yields Eq. 2.4. 

Equation 2,4 for the volumetric reaction rate can be 
written as 


R = n^ np <ov> 


(2.4a) 


where 


<av> 




V a(v) 


q(v;) 


(2.4 b) 


If A and B refer to the same species, f^(v.^) and ^5(^2) 
will be replaced by the same distributions, f(^) and f (^) , 
and Eq. 2,4a will read 


1 9 

R = ^ n <av> , 


(2.4c) 


where n is the total number of particles per unit volume of 
the single species under consideration. The appearance of the 
factor 1/2 in Eq. 2.4c is easily explained on the basis of 
every interaction having been counted twice when considering 
a single species of total number density n [l,6]. Alternatively, 
one can first take = n/2, and consider the interaction 

between these tv;o halves leading to R^ = (l/4) n"^ <crv>. One 
can then add to it successive contributions from within each 
half and so on, finally leading to Eq, 2,4c, 
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2.2 ASSUMPTIONS AND THEIR VALIDITY FOR FUSION REACTIONS 

In deriving Eq. 2,4 for tlie volumetric reaction rate, 

R, various assumptions have heen implicitly made. Some of 
the (more obvious) assumptions, and their validity for fusion 
reaction rate calculations, are briefly discussed belov/, 

(1) Eq. 2.4 assumes that the particles of species A 
and B can be considered as classical particles, i.e, as loca- 
lized wave packets v/hose extensions are small compared to the 
average interparticle distance. This will be realized if the 
average de Broglie wavelength, is much less than the average 
interpjarticle distance, n L12J s 

rJ^'^ << 1 , (2.5) 


X 

d 


iL. 

mv 


(2.58) 


where n is the particle density, h Planck’s constant, m the 

mass of the particles (ions in a fusion reaction) , and v is an 

average speed, say, the rms speed of the particles. If the 

condition 2.5 is not satisfied, the linear dependence of R on 

n^ and n^ (separately) in Eq. 2,4 will breakdown, even if 

suitable definitions of other variables can be arrived at. For 

a fusion reaction at say 10 keV, it is easily seen that X^^ 

10 m which is sufficiently less than the interparticle 

b1 —0 

distance even at n = 10 n , a t 3 rpical particle density for 
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inertial fusion. (Some of the constants needed for these 
estimations are given in Appendix C ) . 

(2) The range of the forces leading to the reaction 
under consideration in Eq. 2.4 should he sufficiently less 
than the interparticle distances. The best way to check it 
is : 


„1/2 „l/3 


( 2 . 6 ) 


—28 2 

For fusion reactions, the highest value of a is 5x10 m , 
and the condition 2.6 is seen to be satisfied even in the 
worst case of inertial fusion. It may be mentioned that if 
condition 2.6 is not satisfied, again the linear dependence 
of R on n^ and n^ will breakdovm, i.e, some other term(s) in 
the rhs of Eq, 2,4 will also have some dependence on n^ and 

ng. 


(3) For Eq, 2.4 to be valid, the velocity of a 
particle should not be correlated to its position with respect 
to another particle. For this to be realized, there should 
be no long-range interparticle forces. This obviously is not 
valid for plasmas. However, if the average potential energy, 
<PE>, in the system per unit volume is sufficiently less than 
the average kinetic enei'gy, <iiE>, per unit volume, the effect 
of correlations on the reaction rate can be ignored [6,4], 

This is equivalent to the condition that the Debye length, 
be large compared to the interparticle distance [2,3j i 
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^2.7) 

where 

T 1/2 

^ ^ » (2.7a) 

n e ■ 

= permittivity of free space, 

T = temperatux'e, in J,and 

e = electronic charge. 

It is seen that for inertial fusion (n'—IO^' m~^, T 10 keV) , 
1 /7 

Xq n ' 1. Hence the condition 2.7 is not satisfied, and 

correlation corrections may be required. However, all the 
calculations in this work relaue to magnetic fusion, in parti- 
cular, to velocity distributions appropriate to mirror systems. 
For magnetic fusion, the condition 2,7 can be easily seen to 
be satisfied, 

(4) It should be mentioned that Eo^. 2,4 gives only 
an expected reaction rates. For magnetic fusion smallest 

— Q ^ 

volme aV and time 5t of interest are of the order of 10 m 

and 10 s, respectively. For inertial fusion, the same may 

be 10 m^ and 10"^ s. At T = 10 keV and densities 

21 — d d1 — d 

of the order of 10 m and 10 m respectively, it can be 

R 

seen that RAV5t • — i 10 . Hence the expected rate ’,;ill be close 
to the actual rate, fluctuations in RAVdt being of the order of 

J RavS [6] . 
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(s) Relativistic corrections are completely negli- 
gible, average ion velocities in fusion plasmas being of the 
6 -1 

order of 10 ms . The effect of external forces on fusion 
cross sections is ignored (Hov.-ever, see the next rem.ark) , 

(6) Polarization effects are ignored, i.e,, in the 
measurement or calculation of o(v) , all allowable relative 
spin orientations of the interacting species are considered 
to be equally probable. ( Oppositely polarized ions can enhance 
fusion cross sections b^^ as much as a factor of tvra [llj. The 
external forces, for example, magnetic fields, do effect the 
state of polarization, generally acting as a depolarizing 
mechanism, and thus destroying the reaction enhancing effect 
of spin polarization.) If the state of polarization is clearly 
known, an appropriately weighted 0(v) can be defined and Eq, 
2=4 for reaction rate calculations can still be used. 

2,3 <av> FOR MMLLIAN DISTRIBUTIONS 


2.3.1 Three Dim.ensional Maxwellian Distribution 


The three dimensional normalized Maxwellian distribution 
of velocities for, say, a species A, whose (peculiar) velocities 
are denoted by, say, v^ , is given by; 


^ m. 3/2 

(^) = ( -^ ) 


2T 


27rT 


f 


( 2 . 8 ) 
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from which it follows that the tliree dimensional Maxwellian 
distribution of speeds is 

As usual, we use the same symbol for the velocity and speed 
distributions, even though the two functions are different. 

In Eqs. 2.6 and 2.8a, m^ is the mass of the particles (in Kg.) 
and T is the thermodynamic temperature (in J). Normalization 
is easily verified by integrating Ea^, 2,8a w.r.t. v^ from 

9 

v^ = o to » . Alternatively, in cartesian coordinates, v^ 

2 2 2 

in Eq, 2.8 can be replaced by v + v + v and integration 

x1 y1 z1 

performed w.r.t. v , v and v from - oo to i- « . in both 

IK. i y ‘ z I 

cases, the result obtained is unity. It should also be pointed 
out that the distribution in Eq, 2,8a is for all solid angles. 

If it is written on the basis of per steradian, the factor Am 
will not appear in the distribution function, but will reappear 
during integration over the entire solid angle. 

In order to calculate <av> (see Eq, 2,4b) for two 
species A and B, both of v/hich are at the same thermodynamic 
temperature T, the relative speed distribution between the two 
species can be -I'jritten down using the basic laws of statistical 
mechanics [6j, i.e, the same laws which yield Maxwellian distri- 
bution of individual species to begin with. This is obtained 
by replacing m in Eq. 2,8a by the reduced mass, 4 : 


L-n ( 

27TT ^ 


3/2 





(2.8a) 
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(2.9) 



(v) 


H 3/2 

4n: ( — -) e 

2TrT 


2T 2 


( 2 . 10 ) 


where v = jv^ ~ | relative speed as in Eq, 2,2. 

The meaning of is that civ gives the munber 

of particles of species B per unit volurae whose relative 
speeds w.r.t. any given particle of species A are between v 
and V + dv. Alternatively, (v) dv can be interpreted 

reciprocally in a similar way. With this interpretation, it 
is then straightforward to rewrite Eq. 2.4b for this particular 
case of tT,vo species being at the same thermodynamic tempera- 
ture T [lj t 


<0V> 


3D 


= 4k (^) 


3/2 


dv v^ o(v) 


4 ' 

2T 




( 2 . 11 ) 


where 3D denotes 3-d.imensional Maxwellian distribution. 

¥e may mention here that the argument leading to Eq. (2.11) 
may appear somewhat subtle. In Section 2.4, we wall obtain 
expressions for' <av> for more general (non-equilibrium) velocity 
distributions. Eq. 2.11 will then aiDpear as a special case. ¥e 
may also mention here that Eq. 2,11 is the standard equation 
for calculating <av> for thermonuclear reactions [l,5j. 
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2 • 3 . 2 One and Tvjo_ Dimensional Kax\^ ;elli Distributions 


Although one and tv/o dimensional Maxwellian distri- 
butions are e:itreme cases and perhaps not of practical 
interest, have calculated <ov> for these distributions 
particularly to see that these limiting results are approached 
smoothly in computer calculations, as v.'ell as in analytical 
formulations e It is interesting, however, to note that if 
these limiting distributions could be approached in practice, 
it can lead to significant reaction enhancing effect in 
cases of practical importance (see Chapter 3 for details). 


As has been remariied above, these limiting cases will 
appear as special cases of the more general expression for 
<av> to be developed in the next section. Our purpose here 
is to indicate that the argument used in Section 2, 3.1 for 
three dimensional Max’/ze Ilian distribution can also be imme- 
diately extended to obtain <ov> for one or two dimensional 
cases. For convenience, v;e denote the temperature T by 
for 1— dimensional case, and by for 2— dimensional case. 

By denoting the speeds by v^^ and respectively for the 

1-D and 2-D cases, we have the follov/ing speed distributions 
for species A i 




, “a , 




( 2 . 12 ) 


[ 1-dimensional I'lax’A/e Ilian distribution 

of speeds] 
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and 




rn 


- 2ti ( 2 ^^ ) 


2 


V 


1r 


(2.13) 


[ 2 — dimensional Maxwellian distribution ] 

of speeds 

By using the same arguments as in Section 2.3.1, one can write 
respective f^g (v) by replacing in Eqs, 2.12 and 2.13 by 
the reduced mass, h, and finally’’ the iollo\;ing expressions for 
<av> can be obtained s 


<av> 


ID 


2 (JL^) ' 

27iT_ 


/•: 


O 


}1Y 

~ 

dv V a(v) e . (2,14) 


and 


l 


, p, . 

<0v>2d = (“m”) _/ 

r o 


2T^ 


dv o{v) e 


(2.15) 


2.4 TRANSFORMATION TO CM AND RELATFAE ’^/ELOCITY COORDINATES 


In order to calculate <av> for two arbitrary velocity 

’"•'Hbf — 

distributions f^( ) and fg( v^ ), Eq. 2.4b has to be inte- 




grated over six variables - 3 for v^ and each. As it 
stands, none of these integrations can be performed without 
making explicit use of a(.v) , because v 


i ^2 " ^1 


contains 


all the six variables of integration. However, by changing the 
variables of integration from v^ and to V and v (defined 


27 


below) , five of the integrations can be performed as o(v) 
does not depend on them. Integration over v can then be 
performed analytically or computationally, making explicit 
use of a(v) , 


¥e define ; 



“'b 



, and 


( 2 . 16 ) 


V 






(2.17) 


It is immediately seen that V is the center of mass velocity, 
and V is the relative velocity, whose nagnituae v v/e have used 
earlier. By changing the variables of integration to V and v, 
Eq, 2,4b can be written as 


<0V> 


> — 7> 

d-^v d^V V a(v) f.(V,v) fr)(V, v) 




-> 

V 2 ) 




V 


) 


where 


h (v. 


. -5s 

d Cv.^, v^) 


(2.18) 


is the magnitude of the Jacobian 


1’ ’'2^ 


i (7 , t ) 


^ (V, v') 

. It should be mentioned here that after the variables 


are changed, in writing Eq, 2.18, we have continued to denote 
the distribution functions by f^ and fg, even though they are 
different functions. Further, the region o;C integration in 2,18 
should corresponding to the region of integration in the original 
Eq, 2.4a. 
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The Jacobian mentioned above is evaluated in Appendix 
A, and its magnitude is seen to be unity. As a result, 



It is seen that in Eq. 2.19, integration w.r.t. d^V can be 
performed in an^^ convenient coordinates. Hov^ever, for inte- 

3 

gration v.r.t. d v,one of the coordinates must be the magni- 
tude of ^ i.e. Vj so that two of three integrations involved 
can be obtained without using a(v) . Spherical coordinates 

3 

for d V will, therefore, generally be suitable. 


2.5 A T¥0-TEf1PERATURE PSEUDO-I#.XW,LLIAN DISTRIBUTION 

In the next section, we will develop an appropriate 
expression for <av> for a chosen non-l-Iaxwe Ilian distribution 
function of speeds. For this purpose, we choose a two- tempe- 
rature pseudo-Haxwe Ilian distribution function, which in some 
sense, can be used to characterize the velocity distributions 
in a mirror fusion reactor. This is based on having one 
temperature, T^, characteristic of the velocity distributions 
along the magnetic field lines, and another temperature, T^> 
characterizing the velocities perpendicular to the magnetic 
field lines. A normalized velocity distribution function of 
this kind for, say, a species A whose velocities are denoted 
"by , can be written as f ollo^irs i 
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2T^ 2T^ 

e e ^ , (2.20) 

where = component of parallel to the field lines, and 

.. 

v^r component of perpendicular to the field lines. 
In Eq. 2,20 if vie put, 


fAC^) 


1 

m. 2 m. 

(-A.) (-A. 

2nT_ 2nT 


= T, and 

2 2,2 

Vi - » 


¥e immediately get the 3-D Maxwellian distribution of Eq. 2,8. 
If we multiply Eq. 2.20 by 27Tv^^dv^^ and integrate from 
Vir = o to oo , we get 1-D Maxwellian distribution of veloci- 
ties ; 


A 



. ^A'<z 

m„ 2 ^ '^z 

( ^ ) e 
2tiT^ 


( 2 . 21 ) 


It is to be noted that the velocity distribution of Eq, 2,21 
above rightly differs from the corresponding speed distribution 
of Eq, 2,8 by a factor of 2. This is because the velocity 
will range from - oo- to + ««■ , 'v/hile the speed v^^ range 

from o to ■«>, . 

Similarly, if Eq, 2.20 is integrated over (from— 

•to + or ) j 2-D Maxwellian distribution of velocities is obtaineds 
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( ) e 

2nl‘ 
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‘A 




( 2 , 22 ) 


which again differs from the corresponding speed distribution 
of Eq. 2,13 by a factor of 27iv^^,. 

2.6 <ov> B"0R THE PSEUDO-KAXiELLIAK DISTRIBUTION 

In order to calculate <av> for tv/o species A and B, 
with velocity distribution as in Eq, 2,20, we first write a 
similar expression for f 3 (*^ )» taking T^ and T^ for species 
B to be the same as that for species A. (if this is not so, 
the method developed in Section 2.4, which is used below, will 
still be applicable. Only the details will be different and 
more involved.) 

In general, in order to change to CM and- relative 
velocity coordinates, Eqs. 2.16 and 2.17 will be solved for 
six components of v^ and v^ in terms of the six. components 
of V and v, and the resulting expressions will be utilized 
in fj^(^) and f 3 (^) to obtain f^Cv,"^'} and fgCv^’^). Then, 
as pointed out earlier, five of the six integrations in the 
general expression for <av> in Eq. 2,19 can be carried out 
without explicit use of o(v) . 

The same procedure is, of course, applicable for the 

distribution being considered. However, the details are 

simplified for this particular case by noting that Eqs, 2.16 
and 2,17 imply for z-component (see Appendix A) ; 
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+ 



(2.23) 


where fj, is the reduced mass as defined earlier and H = m.+m-n, 

A Jd 

Relations similar to Ec. 2.23 1/111 he valid for other compo- 
nents. As a result we also have, 


"^Ir '^2r 


= 4 V, 


2 


M V 


(2.23a) 


Since only v^^ and v^^ occur in the pseudo-Maxme Ilian distri- 
bution of Eq. 2.20 under consideration here, it is easy to see, 
with the help of Eqs. 2.23 and 2.23a, that the Eq. 2.19 for 
<av> reduces to the following for the present case ; 


<av> 


rz 


(“a 

8n\T^ 



m ^ 

2 T 




V? 


2 T. 


d Vd^v V a(v) e 


d 


2T„ 


2 

2T 


(2.24) 


v-'here we have put a subscript rz on <ov> to indicate that this 
applies to the ti/o-temp erasure distribution of Eq. 2.20. 

As promised earlier, five of the six integrations 
involved in Eq, 2.24 can now be performed. This is done by 
writing d-^^V in circular cylindrical coordinates (Fig, la) ; 


d^V = dV^ dO dV^ , 

3/2/ Z ^ 

and integrating with the following limits ; 
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FIG, 1a s Cylindrical coordinates for d''V-integration in Eq.2,24 
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FIG, 1b s Spherical coordinates for d-" v-integration in Eq. 2.24 
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° o to «» 

^2 ' “00 to 30 , and. 

& s o to 2tt (see Fig. la) 

Tv/o of the integrations related to d.'^v can be performed by 
3 

writing d v in spherical coordinates (FIG. 1b) : 

3 0 

d V = V sin9 d9 d& , 

and replacing v and in Eq. 2.24 by : 

= V cos 9 , 

= V sin 9 , 

with limits ; 


9 ; o to 7T 

0 s o to 271 (see FIG. 1b) 

The resulting expression for is z 


<av> 


rz 


3/2 


2 4 
71 ip ij1 /2 
z r 


a(v) 


2T.^ 


I (v) dv 


o 


(2.25) 


where , 


dv^d^-T^) 

O fp rp 

ZL a. _ X.^ 


I(v) = 


dc 


( 2 . 26 ) 


o 
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[a here is related to the original integration w.r.t. 9 by s 
oc = cos 9 I this integration can not he done analytically, but 
it does not require the knov/ledge of o(v)]. 

If we use (= T) in Eq, 2,25, we immediately 

obtain <ov>^jy as in Eq, 2.11, It is also possible to take 
limits in Eq, 2,25 as o, or as ^ o, and obtain 

<cjv>ij-) and <a-v>^-Q of Sqs, 2,14 and 2,15 respectively. However, 
taking these limits is somewhat involved, and the desired 
results are more directly obtained by starting with Eqs, 2.21 
and 2,22, instead of Eq, 2,20, and repeating the steps out- 
'lined in this section for these special cases. 


2.V BASIS OF C0I®>ARIS0N FOR DIFFEFlEHT <av>* s 
\ 

In order xo compare <av>*s for arbitrarily chosen 
velocity distributions ^ fg ( ), we define 
effective (kinetic) temperatures, and Tg, of species A 
and B, respectively, as follows ; 


3 T 

2 A 


3 

2 




“ m^ < v^ >, and 


1 2 

2 niB < v^ > 


For the two-temperature pseudo-Maxwe Ilian distributions [of 
the form given by Eq, 2,20], with the same and for the 
two species, these effective temperatures are easily seen to 


35 


be the same. Denoting this common effective (kinetic) 
temperature b 3 r Tj^, we find by straightforward integration; 


3 

2 


■■K 




(2.27) 


Since is a measure of the average kinetic energy per 
particle, it can be used as the basis for comparing <av> 
values for different distributions. Thus v;hile varying 
and for the two-ternperature distribution, we compare <ov> 
values for a fixed at a time. 


2,8 TRANSFORMATION TO CONVENIENT (DIMENSIONLESS) PARAMETERS 


To facilitate integration of Eqs. 2.25 and 2,26 in a 
more elegant form, and to facilitate some of the further 
discussions, we define 



v/here is the reduced mass, 
ratio parameter, by ; 


(2.28) 

¥e then define the temperature 



and a dimensionless relative speed, by : 


u 


V 



(2,29) 


( 2 , 30 ) 
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These definitions, Eqs, ,2,27 to 2,30, can be introduced into 
Eqs. 2,25 and 2.26. ¥e first introduce Eqs, 2.27 and 2.29 to 
obtain ; 


-(2/3+l)u 


<av> 


I'Z 5 / 3;T ^T 


2 /,u V I 5 / \ ^ 

(tT-) / v^o(v) I(v,Pj e ^ 


K 


. dv , 


(2.31) 


where 


I(v,l3) = 


I 

r 


(1^-2P) .(1-_|5) q ^2 


6 3 T,. 


o 


. da (2,32) 


The above equations bring out the dependence of on Tj^ 

and 3. If we now introduce Eqs. 2,28 and 2.30 in the last two 
equations, the later can be •\,\'ritten as : 


Q‘ 


<av>. 




2 (23+1)^^^ 


(23+1) 


23 


rz . n’ 3 


.Vj^, / a(^o) I On, 3) e 


. d« , 

(2.33) 


where 


(1+23) (1-3) 


I(“,3) 


1 

/ 

o 


1 23 

e 




n 


da 


(2.34) 


It is easily seen that the case 3 = 1 corresponds to <av>^j^ 
of Eq, 2,11, It is also possible to show, with some involvement, 
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that the limits of Eq. 2.31 3.s 3 and p ^ oo exist and 
correspond, respectively,, to and of Eqs. 2.14 

and 2.15, To avoid repetition, appropriate expressions for 
these special cases are direcxly given in the next section, 
where v;e summarize the results of this chapter v/hich are 
actually used in later calculations. 


2 , 9 SUMMARY 


For a two- temperature (pseudo-ilaxwe Ilian) distribu- 
tion of velocities for species A, 




= (- 


m. 






■’■Ahz 


HI 


(_1_) 

27TT_ 


A lr 
2 T„ 


(2.35) 


and a similar distribution, ig( ^2 ) » species B with same 
T and T , v;e have 

o X 


^'2 (2jS+l) 


3/2 


r 


~ J n' p *^K 7 

^ o 


'(23+1) 

23 


(2.36) 


w^here 


1 i:!±23^Ci-:±) ,,2 ^2 


I(w,3) = I e 


. da , (2.37) 
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^ = ^ It 

'-^2 * 2 -^K 


T + •-■ T 
r ^ 2 » 


3 T. 
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K 




w = 


V 


V = 


^2 


, and 


ju. = 


“a ■*■ 


F or 0 = 1 , 


<ov> 


3D 


r 


3 0)2 


6 ^ ^ 

3 / - I 0)-^ o('a)) e do) , 


for 3=0 


<0V> 


1D 



^ Vj^ I 0 ) a(a)) e 
d 


0)-" 

2 “ 


d ,. 


^ » 


and for 3 = ix> f 


<ov> 


2D 


f 


- 0 )" 


2 v-^. I a(w) e do) 


(2.38) 


(2.39) 


(2.40) 


o 



Chapter 3 


NUMERICAL COMPUTATIONS iiL'D RESULTS 

3.1 FUSION CROSS SECTIONS 

3.1.1 Comparison of Expe riment al an d Theoretical Data from 

Ili^L7j..„an^ 

Tables 3,1 to 3.3 sho\./ experimentally measured values 
^DT» ^DD ^ ^ selected energy points. Here the 

incident particle is deuteron,, and the target particle (triton, 
deuteron and He^ respectively) is considered to be at rest in 
the laboratory system. The deuteron energ\^ shovm is, accordingly, 
the energy of the incident deuteron in the laboratory system. 

It is related to the relative velocity, v, of Section 2.9 by 


1 2 
2 “d ^ 


(3.1) 


where m^^ is the deuteron mass. 

It is seen that there is considerable difference in the 
cross sections from the two sources. For all calculations in 
this v;ork we have used the cross sections from Ref. [l]. The 
values in Ref. appear to have been arrived at by consulting 

different sources of experimental data, and were therefore 
considered to be better suited for our puirpose. The detailed 
values are given in appendix B. It should be noted that the 
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cross section is ■the total oi'' 0 £S section for both the 
reactions shown in Eqs. 1.2a aiid 1,2b. 


On theoretical grounds, the cross section for a nuclear 
fusion reaction can be written as [l]s 



(3.2) 


where = energy of the ineracting particles in center of mass 
system, and 

A and B are constants. 

Eq, 3.2 is particularly valid if \i is \;ell below the Coulomb 
barrier between the interacting particles. Theoretical expre- 
ssions of A and B are available ([lj, [Sj). However, consi- 
dering A and B as consta'nts to be empirically determined, a 
least square fit between Eq. 3.2 and experimental data from 
Ref. [l] for DT reaction is shovrn in Fig, 3.1. For this purpose, 
only seven representative experimental values shown in Fig, 3,1 
were used. The unweighted sqxaare error, ^ , in the cross sec- 
tion values, given by, 

^s = S (a - ( 3 . 3 ) 

was minimized, where o is as calculeted from Eq. 3.2 at the 
representative energy points, and is the corresponding 
experimental value. Here it should be mentioned that if A and 




Least squares 
approximation 

cT(E) = Ae^ 

E 

(A = 3687 
B=632) / 



Measured values) 
[REF.d )] 


100 


1000 


Ed (keV) > 

FIG. 3-1 Least squares approximation for apj 
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B are considered as euipirical constants, W in Eq. 3.2 can 
be replaced oy deutercn energy, E^, in the laboratory system. 

It is seen in Fig. 3.1 that thei'e is considerable 
difference between theoretical and experimental values, even 
if the constants A and 3 are determined empirically by least 
squares approxirnaxion. In this con-noction, we may also remark; 

1. In Ref, [sj , a theoretical expression for o{(^) is 
given as follows ; 

- 

o(v/) = a e V , (3.4) 

¥e considered this expression also, but the least squares 
error was more in comparison v/ith thot obtained from Eq, 3,2, 

2. It appears that the values given in Ref, [1J 
actually correspond to a least square approximation between 
experimental data and a plot of In (E^ a(E.^)) vs 

according to Eq. 3.2. This plot is linear as can be seen from 
Eq. 3.2, by replacing W by E^^. ¥e tried to check it. The 
values do not fall on a straightline . (There is some confusion 
here, probably due to a misprint, as in Ref. [l], it is stated 
that the plot In (E^^ d(E^)) vs E^^ is a straightline in accor- 
dance with Eq, 3. 2 1 this is obviously not true.) 

In any case, the cross section values given in Ref,[l] 
are most extensive, and are in reasci able agreement with limited 
experimental data consulted duxing this work. In Ref , [l], <av> 
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values are also given for three dimensioral Maxwellian 
distribution of velocities. Therefore, for the purposes 
of comparison with some oi the calculations reported in this 
work, it was thought best to use the cross section data from 
Ref. [lj, as summarized in Appendix 3. 

3.1.2 Interpolation betwee n Exaerimental VaJLues 

For integrating Eq. 2.36, a^l energy interval of 1 keV 
was used. Since the values given in Appendix B are at inter- 
vals ranging from 1 to 100 keV, sn interpolation had to be 
used. ¥e considered both linear interpolation and an inter- 
polation based on Eq. 3.2, with A and B determined using the 
two end values under consideration. The final results reported 
in this work arc all based on this later interpolation scheme, 

3.2 SOME REMARKS ON COMPUTER CALCUUTIONS 


An outline of the computer program used for calculating 
<av>^^ in accordance with Eos, 2,36 to 2,40 is sho^m in Fig, 
3,2. All integrations v/ere done using Simpson’s rule [l3js 



dx 


“ Lf(x+h) + f(x)j 


(3.4) 


v;ith an error term of the order of Ir' , This was found to be 
adequate when used with an energy interval of 1 keV. However, 
for integrating Eq, 2,37, especially for low values of 3, very 
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FIG, 3, 2 i Outline of the computer program used to 
evaluate <av> 
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small intervals of c were needed to obtain desired accuracy. 

For most calculations 20 to 50 uniform intervals of a over 
the range o to "1 were adequate. But for small values of 
as many as 1000 intervals of a were needed. Further, it was 
found to be the best to merge the e^cponential in Eq. 2.36 v^ith 
the exponential in Eq, 2,37 for the purposes of integration, 
and to use non-uniform, intervals of a. Relatively large 
intervals in the ran^e o to 0,8 and finer intervals in the 
range 0,8 to 1 gave good and fast results. 

The correctness of the calculations \;as verified by 
the follo’.v’ing consiaerations i 

1, The values of <ov>^,^ given by Eqs. 2,36 and 2,37 for 

|3 — and co were calculated and compared with the limit- 

ing results of Eqs, 2.39 and 2.40 respectively. Calculations 
were improved until accur-ate agreement v;as achieved up to 
P = 10"”^ to B = 10^. Since these were the worst cases, the 
accuracy of the calculations for other values of (3 could 
obviously be assumed. 

2. For P = 1, <av> ^ of Eq, 2.36 should correspond to 

"JO jLi 

<av>^j^ of Eq, 2,38. This was checjced. These values were 
further compared with the <av> values given in Ref. [l] for 
3-dimensional Maxwellian distribution of velocities. Expected 
agreement was found, as will be reportec later. 
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3. If we put o(w) = 1, then Eqs. 2.36 to 2.40 can be 
integrated analytically or through the use of tabulated error 
function for j3 = 1 , o and -oo , The results thus obtained 
agreed with the computed results for these cases. Here we 
may point out that actual limits of integration v;ith respect 
to w in Eqs. 2.36 and 2.33 to 2.40 are not o to » , but 
correspond to deuteron energy values from 7 to 1000 keV, as 
the data available v;as limitea to uhis range (see Appendix B) , 

3.3 RESULTS AIJD DISCI SSIOi'S 

<av> - calculauions \/ere dona for the DT, DB and 

3 

DHe fusion reactions for the .-la xv/e Ilian and non-Hlaxwe Ilian 
velocity distributions, ^/ith a view to see the effect on <ov> 
of a departure of velocity distribution from Ilaxv'/ellian. The 
departures considered were limited to those \mich could be 
described by a tv/o-temperature pseudo-Haxwe Ilian distribution 
as discussed in Section 2,5, The values of <ov> for different 
distributions are then compared for the same value of T^^, the 
kinetic temperature, as discussed in Section 2,7. We summarize 
the results in the follov/ing subsections, 

3,3.1 <ov> FOR l-D, 2-D and 3— D KAXW E LLIA-Il Distribution 

The values of <av> available in literature are for 
3-dimensional Maxwellian distributions. These are listed in 
the last columns of Tables 3.4 to 3.6, taken from Ref, [l]. 



Table 
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In the columns just preceding the last columns, are listed 
the results of our calculations [Eq. 2.33j. It is seen that 
there is general agreement. However, in viev; of the fact 
that the input data (Appendix B) used for all our calculations 
is also from Ref, [l], "better agreement could be expected. 

We do not v;ish “o stress this comparison too much, except to 
note that if the^inxjut data is as in Appendix B, then <av>^j^ 
values should be as calculated here. 

Our main interest here is to see the effect of depar- 
ture from Maxvjellian distribution on <av>. For this purpose 
we list in the same Tables 3,4 to 3.6, the calculated values 
of <ov>^j2i and as calculated in accordance v;ith Eqs. 

2,39 and 2,40. As already noted in G’napter 2, these are 
extreme forms of the wo- temperatur-e,pseudoHMaxv;e Ilian distri- 
bution as ^ — > o or oo respectively. This is hardly a situa- 
tion which would ever exist in practice, but establishes the 
upper bounds of the effect on <ov> of possible departures 
from Maxwellian distribution. 

As we see from Tables 3.4 to 3.6, <av>^T^ and 
can be significantly larger than <ov>^j^ upto a certain tempera- 
ture T|^, This temperature is about 10 keV for DT reaction and 
about 30 ke"V for' DD and DHe reactions. In vie\'/ of the fact 
that these are the temperatures of practical interest in thermo- 
nuclear fusion, we arrive at the conclusion that 1-D and 2-D 
Maxwellian distribution of velocities give enhanced reaction 
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rates for all the three fusicr reactions. How'everj at lighter 
temperatures, the reverse effect is seen. 

These effects can be e^pleinsd as fjllows. For the same 
kinetic temperature the average rt-lotive speed is highest 
in 1-D and lowest in 3-1' Mafv.sllian ois-.Tibutions . At lower 
values of Tp , this contributes positively to <jv> as o increases 
v;ith V. But beyond a certain v.o decreases with v. This is 
apparent fr‘om Append!:; ri only for DT reaction. However, the 
same thing happens for DD and DHe^ reactions if v;e push the 
data beyond 1000 keV in Appendix B. In any case, as far as 
numerical calculations in thi S 1 C 4* Iv are concerned, a is taken 
to be zero beyond 1000 keV. 

Thus at higher T^, the gain in v actually reduces <av>. 

3.3.2 <av> F or a Tv /o-Temperature Pseudo-iiaccv/e Ilian 
Dist r ibut ion 

Figs. 3.3 to 3.7 plot of 2,36 as a function 

of the parameter 3 = T /T„. As can be seen from the figures, 

X kl 

O 

while 3 is varied, the kinetic temperatT.ire , % = ^ '^z'^'^r^ 

is kept fixed. The meaning of this is, that as the radial and 
axial temperatures are varied, the average kinetic energy per 
particle is kept fixed, T^- being a measure for the same. This 
has already been discussed in Section 2,7. 

Figs. 3,3 plots the results for the DT reaction at 5 keV 
and for the DD reaction at 30 keV, the temperatures of maximum 




temperatures T|< = 5 and Tk=30l<eV 



at temperatures T|^ = 20 an 
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interesi: in nuclear fusion for these reactions, respeotively. 
In the vicinity of the 3-D Maxv;ellian distribution, enhance- 
ment in the <iv> for DT reaction at 5 keV is seen to be more 
than 205o 0.5 to 2,0) , The enhancement for the DD reac- 
tion at 30 keV is, hov;ever, only about It is seen rhat 

the reaction rate is enhanced if ^ T^, no matter v;hich is 
more. This can also be confirmed by differentiating <av> 

o 

in Eq. 2.36 partially u.r.t. p. It shous an extremum at 3=1, 
as is also apparent frcrn the computer res'<ilts shown in Figs, 
3.3 to 3.7. 

Figs. 3 Ji- and 3.5 show similar data for DD and DHe"^ reac 
tions, at some chosen values of Again an enhancement in 

<av> value due to deviation from kaxwe Ilian distribution is 
seen. 

Fig, 3.6 gives the results for DT reaction at 10 keV, 
another temperature of practical interest for nuclear fusion. 
An enhancement of the order of 10% in the vicinity of 3=1 is 

3 

noted. In this figure, we also give the results for DKe 
reaction at 100 keV and DD reaction at 200 keV. These tempe- 
ratures can be of practical interest in the next generation 
fusion reactors when attairjment of such iiigh tem.peratures 
becomes feasible. It is to be noted, however, that for these 
cases, deviation from 3— D Maxwellian leads to a reduction in 
<av>. The same is the case with the results shown in Fig. 3.7 
This has already been a>:plained qualitatively in Section 3.3.1 
in connection with the Tables 3.^ to 3,6. 
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3.8 "to j) , 1 1 snoM ^ ov> vs v/itih 3 as a parameizer , 
Since <av> changes fast v/ith a narrov; range of temperature 
is cnosen for each plot, so the dilference in <ov> for diffe- 
rent values of |3 can be displayed with clarity. The tempera- 
tures chofjen are uhe ones of most practical interest, e.g, 

10 keV for the DT reaction in Fig, 3,8 and 30 keV for 
the DD and Dhe^ reactions in Figs, 3.9 and 3.10, It can be 
seen that .Lor S=1 , the <ov> is lowest, and it is larger for 
both ^ = 0,5 and 2,0. 3 — ^ o yields the highest value of <ov> 

followed by the value corresponding to jS — . Finally, in 
Fig. 3.11, a representative plot of <av> is is given where 
the reverse is true, 

A suimaary of the important conclusions of this chapter 
together with some suggestions is presented briefly in the 
next chapter. 








Chapter 4 


SIJI#IARY ALT) CUG-GESTICNS 


4 . 1 SUJ'MARY 

After a brief introduction to the present work in 
Chapter 1, we presented a detailed a'rialysis of the problem 
of calculating fusion reaction retes under non-i''Iaxvye Ilian 
conditions in Chapter 2. 

In Sections 2.1 and 2,2, the basic equations were written 
down and the underlying assumptions and their validity v/as 
discussed. The basic equations were then used in Section 2,3 
to calculate <av> for 1~, 2- and 3-dimensional Maxwellian 
distributions. The expressions for these cases are relatively 
easy to establish. To develop a suitable expression for <av> 
for a general two-tempierati.ire pseudoHVIaxwellian velocity 
distribution (Section 2.5), transformation to CM and relativ-'a 
velocity coordinates was done (Section 2.4), The final 
expression was then obtained in Section 2.6. ¥e further 
discussed the basis of comparison for different values of 
<ov> in Section 2.7. In Section 2,8, suitable dimensionless 
parameters were introduced in the various expressions developed. 

In Chapter 3, vai‘ious details of the nmerical computa- 
tions were discussed. The results were presented and discussed 
in Section 3.3, The main conclusions are given below. 
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4.2 CONCLUSIONS 

The calculations reported in this work indicate that 
the fusion reaction rates for the DT, DD and DHe^ reactions 
are enhanced if the velocity distribution deviates from the 
thermodynamic Maxwellian distribution ^ at temperatures of 
maximum practical interest. For er.ample for the DT reaction, 
the reaction rate may be enlxenced by 10 to 20% in the tempe- 
rature range of 5 to 10 keV. Similar conclusion is valid 
for the DKo-^ reaction for 20 to 40 keV range. Enhancement 
in the DD reaction rate is however seen to be rather small 
(see Section 3.3 for details) . 

¥e should stress, however, that the above conclusions 
are based on considering deviations in the form of a two- 
ternperature pseudo -Maxwellian distrioution. Their validity 
for arbitrary deviations will need to be verified by addition- 
al computations and/or analytical techniques. 

The calculations reported here also show that at higher 
temperatures (more than 20 keV for DT and more than 40 keV 

3 . 

for DD and DHe" reactionc) , the deviation from Maxwellian 
distribution reduces the reaction rates, in contrast to the 
opposite effect at lower temperature mentioned above. 

4.3 SUGGESTIONS 

1 . In this work v;e have considered only two-temperature 
pseudo-Maxwe Ilian distributions. Arbitrary departures from 
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Maxwellian distribnticn need to be considered to draw stronger 
conclusions. 

2. Under certain conditions (see Reil[llj), a strong 
non-Maxvre Ilian tail is superimposed on the Maxwellian distri- 
bution, The effect of tliis on the fusion reaction rates can 
be investigated very much along the lines of this work. 

3. The validity of the present calculations breaks 
down at the very high densities encountered in inertial 
fusion due to correlation effects. This should be incorpora- 
ted in the conclusions, [.see Section 2.2j 

4. As was mentioned in Section 2,2^ the polarization 
effects are known to enliance fusion reaction rates by as much 
as a factor of two under soitable conditions. This should 
also be investigated. 

5. For all the calculations v/e have used the data 
shovm in Appendix B. Thera appears to be some need to improve 
this data in view of the discussions in Section 3.1. Further- 
more, we have used a sharp cut-off at 7 and 1000 keV. This 
range could be enlarged by extending the least squares 
approximation curve. This will have some effect on <cv> at 
very low teroperotures (l to 5 koV) or at very high tempera- 
tures (above, say, 100 or 200 koV) . 
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A1 


APPEl^Iy: A 


Proof of d^Vo - d-^'V d"'v Li.e. froni Pq. 2.4 bj 


From Eqs. 2,16 and 2, 17 have, 
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The partsial derivatives in hq. c1 .2 caa~ be substituted 
from Eq. a1,1 to obtain 
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¥e can solve this matrix: by using the nethcci of step by step 
expansion. Then, v/e can find out 




11 ^ _ 4^) (u2 ,Lt2 2 


= - 2 ''2^’ 


m* -MiV 

>1 13 


'iiA - -n ■:2; 

^ in. ^ ^ A 

“a “7 


ci,.^ :- m-:, . - 

-i) 1 j2 

^ ^ rn. -1' m-^ 

R D 


-1 
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A2 




d^V d^v 


Proof of " -‘"^i ''' 

wher^e i _= or z 


Vu ■'■ "'B'^li 


Q -n 

f'iV^ -:■ u vf 

1 ‘ i 


132.1) 


_S-te£S 


From Eqs. 2.16 and 2,17, v;e have 


M V. 


^'li 


u, ^ 


and 


V. 

1 


^2i ^1i 


From where we can have v. . and v,- . as 

1i -1 


V, 


1i 


nv. - Sq V. 
1 -b 1 


and 


^2i 


liV. + m-j- V. 
1 b 1 

“M ' 


Hence, the left hand side of the above nq. a2.1 can be 
written as. 
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m, - mp^v.)^ -s- (MV. 


T )' >- i 


im , M^V^ V - i - ingM^V^ ■;■ 


= . llfvf (raj,-i-m^) + ’' ^ 


2 ^A 

M 7 ^ + V 

1 i'i X 


M vf + a 


= R.E.S. 


2 2 
Tn„ vt . -i- m-io v^... 
1i B ^1 


- -2 2 
M V . + J.t V . 
1 
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APFEI-iDIX B 


valleys of aT^rri, o..r-. and a USED FOR 
DT ’ ijB ^ 

DHe-' 

CALCULAlIOMo REPORTED TN THIS ¥ORR (FRON REF.[lj) 


i 

S.No. ! Deutercn Sner2:yj 
! 

i keV 

j 

a,,oXl02S 

0 

DHe^ 

i 9 

I 

'1. 

7.00 

0.00010 I 

0.0 

’ • 0.0 ; 

2. ! 

- 1 

8.00 

0.00044 I 

0.0 

0.0 

3. ; 

t \ 

9.00 

0.00099 j'; 

0.0 

0.0 

4. i ! 

10.00 

0.0018 '! 

0.0 

0.0 

!' c, 1 i 

ji • ' ! 

0 

0 

« 

0 

CM 

0.055 

0.00054 

0.0 

6. 

50.00 

0.280 i 

0.00222 

0.0 

7. 1 1 

40.00 

0.710 ! 

0.005 

0.000245 

8. ' ■ 

50.00 

1 .40 i 

0.00849 

0.000924 

9. 

60.00 

2 .20 i 

0.013 

0.00225 

iio. ; 

) { 

70.00 

3.05 i 

0.0175 

0.00473 

i11. I 

80.00 

1 3.85 ' 

0.0222 

0.0082 

ji 2 . i 

90,00 

4.45 i 

0.0270 

0.0135 1 

13. ■ 

100.00 

4,80 ■ 

J 

0.0315 

0.0197 

114. 

200.00 

4.00 1 

0.0640 

0.185 

i15. 1 

! 

300.00 

j 

2.30 : 

0.0848 

0.480 

i' 1 6 . i 

i 

400.00 

1 

1,40 i 

0 

• 

0 

00 

0.724 

■17. 1 

500.00 

1,05 I 

0,120 

0.780 

:i 8 . 1 

600 . 00 

' 0,750 

0.135 

0 . 648 

19. 1 

700.00 

0.580 

0.141 

0.490 

i20. 1 

800 . 00 

i 0,450 

0.152 

0.369 

| 21 . ! 

900.00 

0 

04 

(T\ 

0 

0.17c 1 

0.300 

• 

CM 

CM 

1000.00 

i 0.510 

I 

0.180 ! 

0 .244 
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AL-’i'BMLIX C 


TABLE OF SOFIE CONSTANTS AITO CONVERSION FACTORS 

USED IN THIS WORN 


nij^ - 2.0141 u 

m,j, = 3.01605 u 

m ^ = 3.01603 u 
Ile^ 

1 u = 1.660438 X 10"^'^ K 
1 ke/ = 1.602 X 10“^^ u 


h 

6 

, 626 

X 

-.~y-.li. 

10 

J-S 

N = 

0 

6, 

.022 

X 

10^^ i 

:-mol 

e = 

1. 

.602 

X 

10"^^ 

C 

^0 = 

8. 

.854 

X 

10“'^^ 

9 — ' 
C“N 

C = 

r\ 

,998 

X 

10® ms 

-1 


1 barn = 10 


( 1' H '()) t'oil i.ro ! If'i 


7 «' 


[1J 

[ 2 ] 

[3j 

[4] 

[ 5 ] 
[6j 

[7] 

[ 8 ] 
[9J 

[10] 

[11] 


[ 12 ] 

[ 13 ] 


i I ‘ I ‘ ! 
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